Simulations of stable compact proton beam acceleration from a two-ion-species 

ultrathin foil 
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We report stable laser-driven proton beam acceleration from ultrathin foils consisting of two ion 
species: heavier carbon ions and lighter protons. Multi-dimensional particle- in-cell (PIC) simulations 
show that the radiation pressure leads to very fast and complete spatial separation of the species. 
The laser pulse does not penetrate the carbon ion layer, avoiding the proton Rayleigh- Taylor-like 
(RT) instability. Ultimately, the carbon ions are heated and spread extensively in space. In contrast, 
protons always ride on the front of the carbon ion cloud, forming a compact high quality bunch. We 
introduce a simple three-interface model to interpret the instability suppression in the proton layer. 
The model is backed by simulations of various compound foils such as carbon-deuterium (C-D) and 
carbon-tritium (C-T) foils. The effects of the carbon ions' charge state on proton acceleration are 
also investigated. It is shown that with the decrease of the carbon ion charge state, both the RT-like 
instability and the Coulomb explosion degrade the energy spectrum of the protons. Finally, full 3D 
simulations are performed to demonstrate the robustness of the stable two-ion-species regime. 

PACS numbers: 52.40Nk, 52.35.Mw, 52.57.Jm, 52.65.Rr 



I. INTRODUCTION 

In recent years, laser-plasma based ion accelerators 
have drawn increased attention due to their unique fea- 
tures such as compact size, low beam emittance, high 
accelerating gradient, and high peak energy. Several ac- 
celeration mechanisms have been proposed to produce 
such ion beams, for example, target normal sheath ac- 



celeration (TNSA) p}- 
Coulomb explosion 



shock wave acceleration 
and the break-out afterburner 
(BOA) [l3|. In such schemes with laser intensity less 
than lO'^'^W/cm^ , ions within a few wavelengths of the 
foil surface are accelerated by the charge separation field 
Es at the target rear, and shock wave acceleration or 
hole boring [ll| at the target front. If the laser pulse 
is intense enough to evacuate all the electrons from the 
focal spot via J x B force, the ion acceleration is usually 
accompanied by Coulomb explosion or even the thermal 
expansion (l2| . In these CcLSGS, Si linearly polarized laser is 
usually preferable due to its ability to heat the electrons 
through the oscillating part of its ponderomotive force. 

However, both the beam quality and the energy con- 
version efficiency of these mechanisms are far away from 
the requirements of real applications. The inefficiency 
and lack of tunability probably come from the fact 
that the acceleration is only due to space charge ef- 
fects. Recently, with the rapid development of plasma 
mirror technology [l^, both the laser intensity (in ex- 
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cess of lO^^VF/cm^) and the contrast (^ 10~^^) have 
been increased significantly. One of the most straight- 
forward acceleration mechanisms, radiation pressure ac- 
celeration fl^-[T^ (RPA, also called laser piston or 
light sail [iSl), is being re- visited. Circularly polar- 
ized laser pulses are used to reduce the electron heat- 
ing by directly pushing the electrons inward into the 
target In the fist stage of the RPA regime, the 

plasma foil works like a rectifier to convert the trans- 
verse laser field to a strong longitudinal accelerating field. 
The radiation pressure, P = 2I/c, exerted by the in- 
tense laser pulse is up to several tens of Gbar and ac- 
celerates forward the entire irradiated region of the foil, 
where / = 2ag x 1.37 x 10^^ / {X[iJ,m])^{W/cm'^) is the 
laser intensity, c the light speed, uq the laser dimension- 
less electric field and A the laser wavelength. In principle, 
unlimited ion energy gain is obtainable with a sufficiently 
long pulse. 

Much work has been done to investigate the RPA 
mechanism. A semi-analytical RPA model fl3, [3 was 
proposed which is verified by a series of particle-in-cell 
(PIC) simulations. Typically, a "spiral structure" [T7l.[20| 
in phase space is observed, resulting in a mono-energetic 
spectrum with a high energy conversion efficiency. The 
first RPA experime nt |2l| . performed recently, also ob- 
served phase-stable [Tjl ion acceleration from laser radi- 
ation pressure. However, several issues must be solved 
before the scheme is ready for practical applications. 
One of the most important issues is the fast growth 
of an undesirable transverse instability, similar to the 
Rayleigh- Taylor (RT) instability in inertial confinement 
fusion |il4 . ,22J. Some schemes have been suggested to 
suppress the RT-like instability, for example, using a 
properly tailored laser pulse. However, in most schemes 
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the RT instability inevitably tears the foil into many bub- 
bles and clumps [23:]. Though the ions' energy distribu- 
tion is initially monoenergetic, it degrades to a quasi- 
exponential distribution after several laser cycles. In 
Ref. [131, we proposed a new regime of stable proton 
beam acceleration by using a two-ion-species foil. Two- 
dimensional (2D) PIC simulations have shown that the 
two ion species from an initially homogeneous hydrocar- 
bon target can separate into distinct layers at the begin- 
ning of the interaction. The carbon and proton layers 
continue to share an interface, and are accelerated to- 
gether. The high quality proton beam persists after the 
laser-foil interaction ends. 

In this paper, we provide insight into the influences of 
the foil composition and the laser penetration into the 
foil on the ion acceleration. We find that the scheme us- 
ing two-ion-species shaped foils effectively avoids the foil 
deformation and significantly suppresses the RT-like in- 
stability in the proton layer, if the laser intensity and the 
foil composition are appropriately chosen. In the follow- 
ing, we firstly review the ID RPA model and then discuss 
how to extend it to multi-dimensional cases and how to 
further enhance the stable proton beam acceleration. 



II. RPA MODEL AND ID SIMULATION 
RESULTS 

We begin by reviewing the semi-analytical RPA 
model [14, 25]. When an ultra- intense laser pulse irra- 
diates an ultra-thin over-dense foil, the foil can be re- 
garded as a plasma mirror. Most of the laser energy is 
transferred to the foil with only a small part carried away 
by the reflected laser pulse. Assuming the frequency of 
the incident and refiected laser pulse to be wq and wi, 
respectively, we have uji/ojq = (1 — + /3) ~ 1/47^ 

as a result of the red-shift, where /3 is the foil veloc- 
ity normalized by the light speed c and 7 is the rela- 
tivistic factor. As the foil is accelerated to near light 
speed, the laser energy is efficiently converted to ki- 
netic energy of the foil, with energy conversion efficiency 
?7 = 2/3/(H-/3)~l - 1/472 w 1. For an ultrathin foU, 
we may assume that the accelerated region of the foil re- 
mains intact. The foil motion can thus be described as 
follows: 



dm 

dt 
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where p — miUil is the target areal mass density, m^, 
rii, and / are the ion mass, ion density, and foil thick- 
ness, respectively. El represents the laser electric field. 
The foil dynamics in the two-ion-species case is defined 
by the areal mass density p — (mini + m2n2)l, not the 
structure within the foil, which is quite different from the 
coUisionless shock wave acceleration in Ref. By us- 

ing hyperbolic functions and defining /? = tanh9, we can 
get a very compact form of the parameters 7 = coshO 



and 7/3 = p = sinhO. Then, the Eq. ([T]) can be rewritten 
as 



{coshO + sinhBfd{sinhe) = Adt, A 



Er 



2'KCp 



(2) 



Noting \/l + udu = 1^(1-1-^)2 and considering the initial 
condition /3t=a ^ 0, we obtain 
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-7^ = At + -. 
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(3) 



For p^lor/3<C-^;P~ At. As the target approaches 
the speed of light (/3 — 1), p w (jAt)^ . This indicates 
that the ion energy initially increases at a rate of (It/ p)^, 
but slows to {It/p)3 . Further, we can obtain the exact 
solution of the Eq. ([T]) as below 



P^l-Ut) + x{t)]^ -[;{t)-x{t)]-^, 

- m - 1? 



7- 



(4) 
(5) 



w here ^(t) = l/[l + fe(t)2], >c{t) = h{t)/[l+h{t)^]i , ({t) = 
^(1 + /3)/(l - 13), and h{t) ^3At + 2. Eq. g]) is based 
on the assumption that the foil remains intact, i.e. both 
ion species are accelerated to the same velocity. 

We carry out a series of ID simulations to investigate 
the detailed acceleration process by using the PIC code 
VLPL We let the heavier species be carbon ions 

and the lighter species be protons. In the first case, the 
longitudinal length of the simulation box is x = 60A with 
6 X 10'* cells so that the expected density spike can be 
resolved. For simplicity, we take the laser wavelength 
as A = 1.0/im. Each cell contains about 100 numer- 
ical macro-particles. The target is 0.1 A long, located 
at X — 10 A and composed of fully ionized carbon ions 
and protons with the same number density 71.42nc (the 
critical density ric ^ 1.1 x 10'^^ /cm^ for A = l.O^m), 
which corresponds to an electron density rie = 500nc. A 
circularly polarized laser pulse is incident from the left 
boundary at i = 0. The laser intensity follows a trape- 
zoidal profile (linear growth - plateau - linear decrease) in 
time. The dimensionless laser intensity is ag = 100 and 
the duration is tl = I6T0 (~ 52.8/s, ITq - UTq - ITq). 
Absorbing boundary conditions are applied to both the 
fields and particles. 

Fig.[lja) shows the laser intensity evolution. The wave 
front of the laser pulse arrives at the foil surface at t — 
IOTq. We can see that a part of the incident laser pulse 
is initially reflected by the target 'dt t — ISTq because the 
foil is opaque to the laser. The laser energy is transferred 
to the foil until the laser-foil interaction ends at about 
t = 45To. Fig. [T]Jb) presents the distribution of the ion 
density nc, uh-, and the accelerating field E^- At the 
first stage of the RPA, the electrons are pushed out by the 
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J X B force and a strong charge separation field forms at 
the foil rear. Because of the higher charge-to-mass ratio 
Zh/ttih, protons quickly move to the front of the foil at 
the beginning of the interaction. The separation time can 
be estimated by tsep = ^J'^l'rnH / eE^ — 2.5/s. The ions 
then experience different accelerating fields, as shown by 
the red spike in the field distribution. The accelerating 
field inside the carbon layer is much stronger than in the 
proton layer so that the carbon ions can be accelerated 
with the protons. The acceleration process repeats until 
t = 45To, similar to the "snow-plough" in the electron 
acceleration (27| . 

Our simulation results also agree well with the 
Refs. [U HOl, where a typical "spiral structure" was 
observed in a pure hydrogen foil. The fact that both 
the heads of carbon ions and protons interlace with each 
other in phase space, as shown in Fig.[ljc), demonstrates 
the above assumption on the ion acceleration process. 

The averaged proton energy evolution is shown in 
Fig. [Ud). At t = 45To, the proton energy is as high as 
500MeF, which is slightly higher than the carbon ion en- 
ergy 450MeV/u. Such high energies with a pronounced 
mono-energetic peak are yet unreachable in other accel- 
eration mechanisms. For comparison, we also show the 
theoretical predictions of Eg.© in the figure, which agree 
well with the simulation results. 

We also performed simulations to investigate the influ- 
ence of the foil composition on the final proton energy. 
The laser and foil parameters are the same as in the case 
above, except for the ion density ratio. We keep the 
electron density the same but vary the ion density ratio 
nc ■ nn from 1 : 1 (case 1) to 4 : 1 (case 2) and 1 : 4 
(case 3). The areal mass density p in the three cases is 
Pi '■ P\ '■ Pz — 1-06 : 1.00 : 0.86. The ion energy spectra 
at t = SSTq are shown in Fig. [2l Apparently, both car- 
bon ions and protons show a clear mono-energetic peak 
in spite of the different areal mass density p. Overall, 
the carbon ion energy per nucleon is a little lower than 
the proton's. For the case 1 and 2, they have similar 
ion energy spectra due to the close areal mass density. 
For the case 3, the protons can be accelerated more effi- 
ciently because of the lower areal mass density p3. The 
averaged ion energy evolution is also shown in Fig. [Ud) . 
As expected, a similar curve for the case 1 and 2 is ob- 
served though the proton density significantly decreases 
from 71.42nc to 20nc. This demonstrates that the overall 
foil acceleration depends only weakly on the foil composi- 
tion. While reducing the areal mass density, the protons 
tend to be more efficiently accelerated. We also find that 
the average proton energies for the cases are very close to 
the theoretical predictions as marked by the solid curve 
in the figure, though a deviation is observed at the post- 
interaction stage. 

We should mention that radiation reaction (RR) effects 
are generally recognized to become iinportant when the 
laser intensity is above lO^^W/cm^ [H, However, 
recent researches indicate that, for a circularly polarized 
laser pulse, the RR effects become relevant only when the 



foil is thin enough for the laser pulse to break through 
it [sol . [3l|. In this case, the final proton energy can be 
even increased slightly. However, such effects are not the 
focus of this paper and are neglected. 



III. MULTI-DIMENSIONAL EFFECTS AND 
SHAPED FOIL TARGET (SFT) 

The ID numerical simulations above reveal that the 
RPA dominated regime has the potential to generate ion 
beams with favorable qualities such as high intensity and 
short duration. Mono-energetic ion beams with peak en- 
ergy 200 — 500MeV/u are in demand for medical ther- 
apy of deep-seated abdominal tumors, treatment of brain 
cancers |32:], and proton-driven fast ignition of fusion tar- 
gets. 

However, when the ID model above is extended to 
multi-dimensional cases, several issues arise. First, the 
foil is deformed by the incident laser pulse due to the 
transverse laser profile. This leads to strong electron 
heating which degrades the final ion energy spectrum. 
Second, the ultrathin foil is very susceptible to transverse 
instabilities, such as a RT-like instability [H, [1^. The 
RT-like instability is seeded once the laser-foil interaction 
begins, and develops at the unstable interface within a 
few laser cycles. Gradually, the foil surface becomes cor- 
rugated and pierced by the laser radiation and the entire 
target is torn into many clumps and bubbles |2i|]. The 
final energy spectrum then shows a quasi-exponential de- 
cay with a sharp cut-off energy. Third, the ion accelera- 
tion in the RPA regime is usually accompanied by other 
acceleration mechanisms, such as TNSA, Coulomb explo- 
sion, and even thermal expansion, making the accelera- 
tion process complex and nonlinear. 

In order to avoid foil deformation, we employ a shaped 
foil target (SFT) (33l - [35| to compensate for the trans- 
verse laser profile, as suggested by the Eq.(IT]), i.e., 
i?£ (X p cc I. Considering the usual Gaussian laser 
pulse, for example, the transverse foil thickness ly fol- 
lows ly = 'max[loexp{—y'^ /a'^),li], where Iq is the max- 
imal foil thickness, li the cut-off thickness, and (Jt the 
laser spot size, as shown in Fig. [3] Such a foil geometry 
can be fabricated for experiments by employing polyethy- 
lene foils [2l|. Alternative methods include molecular 
beam epitaxy technique (MBE) fsE] or deposition tech- 
niques for thin films, such as physical vapor deposition 
(PVD) [37l or chemical vapor deposition (CVD) [Hi. By 
directly matching the transverse laser profile and target 
areal mass density, each foil element Ixdly behaves like 
in the ID case and the irradiated spot can be uniformly 
pushed forward. A single ion species SFT irradiated by a 
circularly polarized las er p ulse has been numerically in- 
vestigated in Refs. |33l43a |. where an improved accelera- 
tion structure was shown. However, transverse instability 
of the foil still degrades the initially mono-energetic ion 
beam. In the following, we show 2D simulations of two- 
ion-species foils, i.e., hydrocarbon SFTs, which present 
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difiFerent dynamics from the hydrogen SFT in Ref. (33| . 



IV. 2D SIMULATION RESULTS 

Two basic cases are first investigated with 2D simu- 
lations. In the case A, the simulation box is X x Y = 
50A X 50A, sampled by 10000 x 5000 cehs. Each ceU con- 
tains 100 numerical macro-particles in the plasma region. 
The foil is initially located at x = lOA with parameters 
Iq = O.IA, li = 0.05A, and ctt = 7A. Both species 
have the same particle density nc = nn = 45.7l7ic, 
which corresponds to the electron density rie — S20nc 
3.5 X 10^'^/cm'^). A circularly polarized Gaussian 
laser pulse with the focal size cr^ = 8A is incident from 
the left boundary. The laser duration is = IOTq 
33/s, ITo — STq — ITo). All the other parameters are 
the same as in the ID cases above. 

Fig. m shows the detailed ion acceleration process. As 
expected, the foil deformation is avoided effectively. The 
central part of the foil with the transverse size ~ can 
be accelerated as a whole. We again find that the pro- 
tons move to the front of the heavier carbon ions at the 
very beginning. This can be understood by considering 
the different charge-to-mass ratios of protons and carbon 
ions, and the different spatial location of equilibrium for 
each ion species within their respective effective poten- 
tials in the accelerating frame [25]. Such a separation 
is determined by the balance between the inertial and 
the electrostatic forces, as shown in Fig. [5ja). However, 
in comparison with the ID case we observe additional 
effects such as the RT-like instability and transverse fil- 
amentation which effect the ion acceleration. As a result 
of transverse effects, some carbon ions spread into the 
space trailing the target after t = 22.5To and evolve into 
a cloud instead of a compact bunch as seen in Fig. ^a). 
However, a dense front always trails the protons. The 
sharp front separating the two species is well defined, as 
shown in Fig. [H^a) . 

We observe a clear RT-like instability in the 2D sim- 
ulations. A typical "A structure" [1^ can be recognized 
in both the ion and electron density distributions, as 
shown in Fig. 31 The proton layer is strongly effected, 
which can be understood by considering the relativis- 
tic laser transparency of the foil, which occurs when 
a > 7T{ne/jnc){l/X). In view of the ultralow thickness 
of the foil in this case, the electromagnetic wave can 
partially penetrate the foil and interact with both ion 
species as seen in Fig. ^i) . This would seed the RT in- 
stability at the proton-carbon interface. Gradually, the 
protons also suffer from the RT instability and show an 
obvious "A structure". However, the trailing carbon ion 
front contributes to slowing down the local bunching of 
the protons and hence suppresses the transverse proton- 
RT instability. This is the reason why we still observe a 
semi stable layer of protons in front of the carbon ions 
at a later time, i.e., t = 30To, though the surface of 
the layer is strongly corrugated. Fig. ^c) shows the ion 



energy spectra at t=20, 30, and 40To. Although a mono- 
energetic peak is observed at the very beginning, it be- 
comes broader after t = SOTq. Eventually, the protons 
also evolve into a cloud in space. 

In view of the negative effects of the laser penetration, 
we propose using a moderate intensity laser pulse or a 
higher electron density foil to avoid the RT-like insta- 
bility from effecting the proton layer. We estimate the 
critical laser intensity for this case by 

a ~ TT -,nei = Zcxnc (6) 

Uc A 

For comparison, we perform another simulation to 
show improved proton beam acceleration. 

In case B, the foil is composed of fully ionized carbon 
ions and protons with number density ratio nc ■ nn = 
4:1. The electron density is increased to rie — SOOrtc 
5.5 X lO^'^/cm'^) so that the foil is opaque to the inci- 
dent laser pulse. All the other parameters are the same 
as in the case above. Fig. [B] shows the simulation results. 
In this case, the foil is well maintained for a much longer 
time as compared with the case A. The protons again 
form a distinct layer from the carbon ions and always 
ride on the carbon ion front. The laser pulse does not 
penetrate the foil and the radiation pressure mainly acts 
on the carbon ion layer, as depicted in Fig. (Hl^e-f). As a 
consequence, the proton layer is kept stable and is less ef- 
fected by the evolution of the RT instability in the carbon 
layer. Compared with case A, the sharp front separating 
the two species is much more smooth. Most electrons are 
moving together with the carbon ions, and the spread 
carbon layer act as a "buffer" for proton acceleration. 
This leads to a more stable acceleration structure, which 
is shown by the energy spectra in Fig. ^d). We can see 
that the energy peak is more pronounced and the spread 
(FWHM) is about 20% at t = 30To, which is only half 
of that in the case A. Meanwhile, we find that the cut- 
off energy of the carbon ions always peaks the proton 
spectrum in the both case A and B. 

It is interesting to note that a similar phenomenon was 
also observed in Ref. [s^ \^ , where mass-limited targets 
were studied and the results were interpreted as the di- 
rect Coulomb explosion (DCE). In that regime, electrons 
are overtaking ions as flying compact layers and the re- 
turn current is inhibited when mass-limited targets are 
used. The accelerating field from the Coulomb explo- 
sion of the heavier ions accelerates the lighter protons 
forward. In our case, however, the laser radiation pres- 
sure is dominating the proton acceleration because the 
electrons always accompany both the carbon ions and 
protons. The bulk electrons enable a good screening of 
the electrostatic field and thus effectively prevent the ions 
from the Coulomb explosion. The foil is accelerated for- 
ward as a compact dense quasi-neutral plasma until the 
carbon ions evolve into a cloud in space. However, we find 
that the Coulomb explosion of the proton layer becomes 
significant later. The peak proton energy and the cut- 
off energy of the carbon ions continue to increase slowly 
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after the laser-foil interaction ends. After the interac- 
tion, i.e., t — 35To, the carbon ions evolve into a cloud 
in space and the electrons are extensively heated, and 
expand with the carbon ions. Finally, the Debye length 
of the electrons may become larger than the proton layer 
thickness. Assuming the thickness of the proton layer be- 
ing Ih.o initially and In^t at a specific time point, we have 
TTcrj^riH.olHfl — ''^cr^'nH,tlH,t according to particle number 
conservation. We can then get a condition for the onset 
of the Coulomb explosion in this two-ion-species regime 
as follows 



T,^t[MeV]>8x{ 



.A 



(7) 



where Tg.f and ng.t are the electron temperate and den- 
sity in the proton layer at the time t, respectively. At 
the post-interaction stage, the electron heating is signifi- 
cant so that the electron temperature is much larger than 
the value by Eq. ([7|), and the Coulomb explosion occurs. 
For example, at t = 35Tq in the case B, the proton layer 
thickness is about Ih = 1/xm and the averaged Debye 
length of the electrons in the proton layer is d = 1.6/im. 
In this situation, the effects of the Coulomb explosion 
should be taken into account (4ll . |42| . Coulomb explo- 
sion is undesirable because it broadens the final proton 
energy spectra, though the ion energy can be increased 
slightly 

It is also necessary to note a recent paper by Qiao, et 
al., [43] , where they also observed the ion species sepa- 
ration by the laser radiation pressure in the "leaky light 
sail" regime. However, we did not see such a wide spatial 
gap between the heavier carbon ions and the light pro- 
tons. Both the simulations by VLPL and VORPAL Q 
show almost the same results as elucidated above. 



V. THREE-INTERFACE MODEL AND 
DISCUSSIONS 

The stability of the proton acceleration in the 2D cases 
above can be attributed to two crucial effects. First, the 
protons advance ahead of the carbon ions and form a 
thin layer. Such a separation of the ion species can be 
understood within the ID formalism developed in Sec. (II) 
and the ID simulations above. Second, a much thicker 
and heavier ion layer trailing behind the proton layer 
acts as an effective "buffer" to prevent short- wavelength 
perturbations from feeding through into the thin pro- 
ton shell. As was shown earlier [24], in the absence of 
such a buffer (e.g., when carbon ions are either absent or 
spatially separated from the protons), the entire proton 
layer is affected by the RT-like instability. We introduce 
a simple three-interface model, as shown in Fig. [71 to ex- 
plain the stabilization of the proton beam acceleration. 
In our case, both species have two interfaces: one with 
vacuum and one with the other species. The only un- 
stable interface is the carbon-vacuum boundary, where 



the laser radiation pressure directly acts on the carbon 
plasma. Below we have basic single-fluid hydrodynamic 
equations describing the RT instability. Detailed deriva- 
tion of these equations will be presented elsewhere. The 
following assumptions are used: (a) electron temperature 
is much higher that the ion temperature: Tg S> T; (b) 
plasma ions are fully neutralized by the plasma electrons; 
(c) the entire plasma is accelerated with the constant ac- 
celeration Under these assumptions, the single-fluid 
momentum equation in the accelerated in the z— direc- 
tion reference frame can be written as: 



Pi 



dv 



V ■ Vw) — — Vp — pg. 



(8) 



where p, v, and p are the mass density, velocity, and ion 
pressure, respectively. Linearizing this equation around 
the plasma equilibrium, and assuming an isothermal elec- 
tron response such that dp — {Zi/mi) {5p)T(., we find 
that the perturbed pressure 6p satisfies 



dv —J. ^ _ 
p— = - Vdp - dpg, 



(9) 



which can be even further simplified by assuming incom- 
pressibility of the plasma in response to the RT-like in- 
stability: 



k%r^Sp 



(10) 



where kuT is the perturbation wave-number in the di- 
rection perpendicular to the acceleration, i.e. paral- 
lel to the interface. Noting that 6p is discontinuous 
across the unperturbed boundary, we obtain a solution 
Sp — Ae~^'^'^^' -t- Be^"'^^ away from interfaces, with A 
and B being the amplitude coefficients of the perturba- 
tion. These solutions have to be matched across every 
interface to satisfy the following continuity conditions: 
(i) Vz is continuous (consequence of the single-fluid de- 
scription); (ii) Sp — i^poVz is continuous (consequence of 
the pressure balance across the perturbed boundary) . We 
can finally derive from this simple model that the ampli- 
tude of the RT instability is exponentially decaying away 
from the unstable interface: 



V. 



H/C 



C /vac 



(11) 



where k is the thickness of the "buffer" ion layer, and 
v^^'^ iv^J^"''^') are the velocity perturbations at the pro- 
ton/carbon (carbon/ vacuum) interfaces. In the case B, 
we know that Ic is a few times larger than In. In this sit- 
uation, the long-wavelength perturbation in the carbon 
layer should take much more time to grow. Finally, the 
perturbation is exponentially attenuated before reaching 
the proton layer. We can interpret the carbon ion cloud 
as a "buffer" or "cushion" for the proton- RT instabil- 
ity. The carbon ion front always trails the protons and 
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helps to confine the spatial development of the proton- 
RT instability, i.e., local bunching of the proton beam ("A 
structure"). Thus the RT-like instability in the proton 
layer is significantly suppressed. 



Table I: Case List 



Label 


Composition nc'-nn 


Tie 


Case A 


C^+, H+ 


1:1 


320nc 


Case B 


C^+, H+ 


4:1 


500nc 


Case C 


C«+, D+ 


4:1 


500nc 


Case D 


T+ 


4:1 


SOOric 


Case E 


H+ 




500nc 


Case F 


C5+, H+ 


4.8:1 




Case G 


(7^+, H+ 


8:1 


500nc 



It is helpful to consider the pro blem from the view 
of the classic RT instability |45l l46l| which occurs when 
a light fluid is accelerated into a heavier fluid. In the 
simulations, the proton layer is much lighter than the 
subsequent carbon layer, and thus proton layer can keep 
stable. We also perform simulations where we employ 
carbon-deuterium (C-D foils, case C, Zo/mE, = 1/2) and 
carbon-tritium (C-T foils, case D, Zt jmT = 1/3) instead 
of hydrocarbon (C-H foils, case B, Zn/mH — 1/1). AH 
other parameters are the same as in the case B. The sim- 
ulation results are depicted in Fig. |8l As discussed in 
Sec. (II), the separation of ion species depends on the ion 
charge-to-mass ratio. For the C-D foil, both the carbon 
ion and the deuterium have the same Zi/rrii and there- 
fore do not separate from each other (see Fig.|8Kb)). The 
laser pulse pushes both of them forward together. In this 
case, it is expected that the RT instability shall deteri- 
orate both ion energy spectra (see Fig. IHd)). For the 
C-T foil, we do observe ion separation initially because 
the carbon ions have a larger Zi/nii and hence they move 
to the front of the tritium. However, the tritium is ini- 
tially compressed into a very thin layer, which leads to 
fast growth of the perturbation at the vacuum-tritium in- 
terface. The instability soon reaches the tritium-carbon 
interface and pollutes the front carbon ion beams. Obvi- 
ously, such a structure is unstable, and we do not observe 
a mono-energetic beam as seen in Fig. |8l[d). 

We also investigate effects of the charge state of the 
carbons on the stability of the proton acceleration. In the 
simulations, we keep the proton density same but vary 
the charge on carbon ions. (case F, Z/m = 1/2.4) 

and (case G, Z/m — 1/4) are taken into account. 
For comparison, a pure hydrogen SFT (case E) is also 
considered. All the other parameters are the same as 
that of the case B. The simulation results and the proton 
energy spectra are presented in Fig. [51 

It is interesting that with the decrease of the carbon 
ion charge state, the gap between the carbon ions and 
protons grows, which leads to poorer monochromaticity 
of the proton energy spectrum. It implies that the RT 
instability is much more severe in the case of the lower 
charge state of the carbon ions. This effect is beyond the 



scope of the current three-interface model and will be 
studied further. Besides, we believe that the Coulomb 
explosion of the proton layer also plays a role in case G 
because most electrons are with the carbon ions, far away 
from the protons. The wide gap prevents the electrons 
from neutralizing the proton layer. Finally, the electron 
density becomes very low in the proton layer. Accord- 
ing to Eq. ([7]), the electron temperature for the onset 
of the Coulomb explosion is reachable. For example, at 
t = 20To, the proton layer thickness, Ih = 0.3/im, is 
smaller than the averaged Debye length of the electrons 
in the proton layer, Ad = 0.4/j,m. As shown in Fig. IH^c), 
we observe a clear dipped accelerating field with a nega- 
tive component in the proton layer, which is a clear signa- 
ture of Coulomb explosion. Gradually, the protons self- 
expand in space and the energy spectrum broadens. This 
is different from the case B where the Coulomb explosion 
might occur after the interaction ends. Additional simu- 
lations show that when the charge state of the carbon ions 
is lower than 3, both the RT instability and the Coulomb 
explosion of the proton layer become violent enough to 
destroy the mono-energetic proton beam. However, in 
the relativistic regime, the laser pulse is powerful enough 
to fully ionize the carbon atoms. We therefore believe it 
is possible to observe stable proton acceleration in the in 
experiments. 

Besides, we believe that a smaller transverse size of the 
foil would benefit the stabilization of the proton acceler- 
ation in this two-ion-species regime. In the simulations, 
we apply a periodic boundary condition in the transverse 
direction and use the transverse size of the simulation box 
Y = 50A. When we decrease this size toY = 25A, we ob- 
serve considerable suppression of the penetration of the 
perturbation into the target. 



VI. 3D SIMULATION RESULTS 

Finally, 3D simulations are performed to check the ro- 
bustness of the stable acceleration mechanism. To save 
computational time, the 3D simulation box is 40A x 25A x 
25A, sampled by 4000 x 200 x 200 ceUs. Each cell con- 
tains 27 particles initially. The ultrathin SFT is com- 
posed of fully ionized carbon ions and protons with the 
density ratio nc ■ n-n =4:1. Periodic boundary con- 
ditions for particles and absorbing boundaries for fields 
are applied. The total laser energy is ~ 500J and the 
duration is 33/s. All the other parameters are the same 
as in the case B above except ax = 5A to match the 
laser focal spot ctl = 6A. Fig.fTUTa.b) illustrates the den- 
sity distributions of the protons and the carbon ions at 
t = AOTo . We can see a clear compact proton bunch with 
a few nano- Coulomb riding on the front of the carbon 
ions, which agrees well with the 2D simulation results. 
The carbon ions spread extensively in space. This should 
be attributed to multi-dimensional effects, i.e., the fast 
growth of the RT instability in the carbon ion layer, as 
discussed in Sec.(V). 
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In Fig. fTUT c). we present the ion energy spectra at 
t ~ 40To. As expected, we observe a clear mono-energetic 
peak for the protons. The peak is well maintained for a 
very long duration even after the laser-foil interaction 
ends. On the contrary, the spread carbon ions show 
a quasi-exponentially decaying in the energy spectrum. 
The cut-off energy of the carbon ions is almost the same 
as the peak energy of the protons, as observed in the 2D 
cases. It indicates that the carbon ion front indeed trails 
the proton layer during the acceleration. We also find 
that there is another small energy peak for protons in 
the energy spectrum. By analyzing the ion energy dis- 
tributions in space, we find that the protons with energy 
~ 280MeV are situated in the vicinity of the carbon ion 
front. We believe it is due to the incomplete separation 
of the protons from the heavier carbon ions in the 3D 
case. Fig. Iior d) presents the proton peak energy evo- 
lution both in the 2D and 3D simulations. Overall, 3D 
simulation results fit well with the ID model marked by 
the red curve and the 2D simulations by the green curve. 
This demonstrates the domination of the RPA mecha- 
nism. The yielded proton beam with the peak energy 
^ AOOMeV/u may have diverse potential applications in 
the future, such as in medical therapy of deep seated 
tumors and in the development of future compact ion 
accelerators. 



VII. CONCLUSIONS 

In conclusion, we investigate the detailed ion accel- 
eration from an ultra-thin hydrocarbon foil by use of 
multi-dimensional PIC simulations. A stable compact 
proton beam acceleration is observed, for the first time, 
in the 3D geometry. This should be attributed to two 
effects: ion species separation and heavier ion spread- 



ing in space. The laser pulse does not penetrate the foil 
and the radiation pressure mainly acts on the carbon ion 
layer. The carbon ions act to buffer the compact proton 
layer from the RT-like instability. The proposed three- 
interface model well describes the simulation results and 
is further supported by simulations of various compound 
foils such as C-D, C-T, and pure hydrogen foils. It is also 
found that with the decrease of the carbon ions charge 
state, both the RT instability and the Coulomb explo- 
sion become increasingly violent and tend to degrade the 
mono-energetic proton beam. Finally, the robustness of 
the stable two-ion-species regime is checked by the full 
3D simulations. 

With the development of the nano-technology (47| . 
polyethylene or mylar foils could be available soon in ex- 
periments. Compared with the normally used two-layer 
foils, they are much easier to fabricate. Benefiting from 
new state-of-the-art lasers such as HiPER and ELI (isj . 
we believe that the stable acceleration mechanism de- 
scribed above will be experimentally demonstrated soon 
and has the potential for applications in science and 
medicine. 
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Figure Captions 

Fig. 1. (Color online) (a) Laser intensity evolu- 
tion. Here, I = {Ey + Bl)/2 is normalized to /q = 
2.74 X lO^^W/cvn?. The laser-foil interaction ends at 
t = 45To. (b) Carbon ion and proton density distri- 
butions at t = 25To and 45To. The red color shows 
the accelerating field which is normalized to Eq = 
ruecuj/e = 3.2 x (c) Carbon ion and proton 

phase space distributions at t = 25To and 45To. (d) Pro- 
ton energy evolution from ID PIC simulations and ID 
RPA model. 

Fig. 2. (Color online) Energy spectra of (a) carbon 
ions and (b) protons a.t t = 35To. The ion density ratio 
nc ■ uh in all three cases is 1:1 (case 1), 4:1 (case 2), and 
1:4 (case 3), respectively. A clear mono-energetic peak 
for both carbon ions and protons is observed in all three 
cases. 

Fig. 3. (Color online) Schematic of the two-ion-spices 
shaped foil target (SFT). A Gaussian laser pulse is in- 
cident from the left boundary. The profile of the foil is 
defined by three parameters: Iq, li, and ax- The foil is 
composed of mixed carbon ions and protons with various 
particle density ratios nc ' nn- 

Fig. 4. (Color online) Density distributions of protons 
and carbon ions in the case A at (a) t = IT.STq and (b) 
22.5To. RT instability is observed initially for both car- 
bon ions and protons. The corresponding electron den- 
sity distribution and the laser electric field Ez are shown 
in frames (c-d) and (e-f), respectively. Here, the density 
is normalized to ric and the field to Eq. 

Fig. 5. (Color online) (Color online) Particle density 
distributions and the electron energy distribution as well 
as the accelerating fields E^ on the laser axis a.tt = 20Tq 
in the (a) case A and (b) case B. Energy spectrum evo- 
lutions of the protons and carbon ions at t = 20, 30, 40ro 
in the (c) case A and (d) case B. For the case B, the 
compact proton bunch can be accelerated steadily and 
the beam quality is well maintained as time goes on. 

Fig. 6. (Color online) Density distributions of protons 
and carbon ions in the case B at (a) t = 17.5To and (b) 
22.5To. The proton-RT instability is significantly sup- 
pressed as compared with the case A. The corresponding 
electron density distribution and the laser electric field 
Ez are shown in frames (c-d) and (e-f), respectively. 

Fig. 7. (Color online) Schematic of the three-interface 
model. The top red marks the protons and the bottom 
blue shows the carbon ions. The laser pulse is irradiating 
on the foil from the bottom. Three interfaces can be 
recognized in the figure: carbon-vacuum, carbon-proton, 
and proton- vacuum. 

Fig. 8. (Color online) Ion density distributions in the 
(a) case B (C-H foil), (b) case C (C-D foil), and (c) case 
D (C-T foil) at t = SOTq. In the C-D foil, both species do 
not separate from each other, (d) Ion energy spectra in 
the case B, C, and D at t = 40To. Here, C(B), C(C), and 
C(D) denote the carbon ions in the case B, C, and D, 
respectively. A stable acceleration structure exists only 



in the case B. 

Fig. 9. (Color online) Density distributions of protons 
and carbon ions ((a) case E,(b) C^"*", case F, and (c) 
C^+, case G.) at t = 35To. (d) Proton energy spectra in 
the case B, E, F, and G at t = 35To. 

Fig. 10. (Color online) Density contours of (a) protons 
and (b) carbon ions in 3D simulations at t = 40To. A 
clear proton bunch is formed before the carbon ion front. 
The energy spectra of protons and carbon ions are shown 
in the frame (c) . The proton peak energy evolution in the 
2D case and 3D case is shown in frame (d). For compar- 
ison, the red color marks the theoretical predictions of 
the ID RPA model. 
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Figure 1: (Color online) (a) Laser intensity evolution. Here, 
/ = {E^ + B^^)/2 is normalized to 7o = 2. 74x 10^® W/cm^. The 
laser- foil interaction ends at t = ibTo. (b) Carbon ion and 
proton density distributions a.t t — 25To and 45ro. The red 
color shows the accelerating field Ex which is normalized to 
Eo — nieCuj/e — 3.2 x lO^^V/m. (c) Carbon ion and proton 
phase space distributions a.t t = 25To and 45ro. (d) Pro- 
ton energy evolution from ID PIC simulations and ID RPA 
model. 




Figure 2: (Color online) Energy spectra of (a) carbon ions 
and (b) protons a.t t = SSTq. The ion density ratio nc : riH 
in all three cases is 1:1 (case 1), 4:1 (case 2), and 1:4 (case 
3), respectively. A clear mono-energetic peak for both carbon 
ions and protons is observed in all three cases. 
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Figure 3: (Color online) Schematic of the two-ion-spices 
shaped foil target (SFT). A Gaussian laser pulse is incident 
from the left boundary. The profile of the foil is defined by 
three parameters: Iq, h, and or. The foil is composed of 
mixed carbon ions and protons with various particle density 
ratios nc : nn. 
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Figure 4: (Color online) Density distributions of protons and 
carbon ions in the case A at (a) t = IT.STb and (b) 22.5To. 
RT instability is observed initially for both carbon ions and 
protons. The corresponding electron density distribution and 
the laser electric field Ez are shown in frames (c-d) and (e-f), 
respectively. Here, the density is normalized to ric and the 
field to Eq. 



12 




rj.5 I3.S I4.I t4J 14.7 15.(1 

X/X 



^ ' Case A 


t-20T(i« 

— t=30Tii, H' 
- t=40T||, H" 

— t=40Tii, C" 








12.6 12.9 13.2 13.5 13.: 

X/?.. 



2U0 400 600 800 1000 

E/MeV 



Figure 5: (Color online) Particle density distributions and the 
electron energy distribution as well as the accelerating fields 
Ex on the laser axis at t = 2QTo in the (a) case A and (b) case 
B. Energy spectrum evolutions of the protons and carbon ions 
at t = 20, 30, 40ro in the (c) case A and (d) case B. For the 
case B, the compact proton bunch can be accelerated steadily 
and the beam quality is well maintained as time goes on. 
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Figure 6: (Color online) Density distributions of protons and 
carbon ions in the case B at (a) t = 17.5To and (b) 22.5Tb. 
The proton- RT instability is significantly suppressed as com- 
pared with the case A. The corresponding electron density 
distribution and the laser electric field Ez are shown in frames 
(c-d) and (e-f), respectively. 



13 



Vacuum 



stable 




Vacuum 



Figure 7: (Color online) Schematic of the three-interface 
model. The top red marks the protons and the bottom blue 
shows the carbon ions. The laser pulse is irradiating on the foil 
from the bottom. Three interfaces can be recognized in the 
figure: carbon-vacuum, carbon-proton, and proton-vacuum. 
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Figure 8: (Color online) Ion density distributions in the (a) 
case B (C-H foil), (b) case C (C-D foil), and (c) case D (C- 
T foil) at t = SOTq. In the C-D foil, both species do not 
separate from each other, (d) Ion energy spectra in the case 
B, C, and D at t = 40To. Here, C(B), C(C), and C(D) denote 
the carbon ions in the case B, C, and D, respectively. A stable 
acceleration structure exists only in the case B. 
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Figure 9: (Color online) Density distributions of protons and 
carbon ions ((a) case E,(b) C^"*", case F, and (c) C^+, 

case G.) at t — SSTq. (d) Proton energy spectra in the case 
B, E, F, and G at t = 35To. 




Figure 10: (Color online) Density contours of (a) protons and 

(b) carbon ions in 3D simulations at t = 40To. A clear proton 
bunch is formed before the carbon ion front. The energy 
spectra of protons and carbon ions are shown in the frame 

(c) . The proton peak energy evolution in the 2D case and 
3D case is shown in frame (d). For comparison, the red color 
marks the theoretical predictions of the ID RPA model. 



